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fractions  with  those  obtained  from  chemical  analysis  was  within  15  percent  for 
binary  component  cryodeposits ,  but  much  worse  for  multicomponent  cryodeposits. 
During  the  search  for  an  explanation  of  the  discrepancies,  it  was  found  that 
the  accuracies  of  the  densities  of  the  cryodeposits  were  as  crucial  as  those 
of  the  band  strengths,  and  that  the  above  assumption  roust  be  modified,  in  a  yet 
unknown  manner,  for  materials  where  there  is  a  large  degree  of  hydrogen  bonding. 
Also,  it  was  noted  that  in  mole  fraction  determinations  the  use  of  absorption 
index  values  alone  is  limited  to  materials  having  at  most  one  IR-inactive  molec¬ 
ular  species.  However,  the  identity  of  IR- inactive  species,  which  must  be 
known  in  the  calculation,  must  be  available  from  techniques  other  than  IR 
spectroscopy.  A  method  using  molar  ref ractlvities  for  analyzing  materials  with 
two  or  more  IR- inactive  species  is  outlined. 

Kramers-Kronig  and  subtractive  Kramers-Kronig  computer  programs  were 
written  for  the  analysis  of  non-normal  reflectance' data  from  a  single  interface, 
for  the  calculation  of  absorption  index  values  from  an  n-spectrum,  and  for  the 
finding  of  the  optical  constants  of  a  single  film  thickness  on  a  thick  substrate 
using  transmittance  measurements.  A  Kramers-Kronig  program  for  obtaining  the 
optical  constants  of  a  single  film  thickness  also  seems  possible.  The  need  for 
more  optical  data  of  substrate  materials  at  cryogenic  temperatures  is  stressed,  i 
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PREFACE 

The  research  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC).  The  analysis  was  performed  by 
Westminster  College,  Fulton,  Missouri,  under  subcontracts  No.  78-22-VKF  and  79-18-VKF 
for  ARO,  Inc.  (a  Sverdrup  Corporation  Company),  operating  contractor  of  AEDC,  AFSC, 
Arnold  Air  Force  Station,  Tennessee.  The  work  covered  the  period  from  October  1 ,  1977  to 
September  1,  1979  and  was  done  under  ARO  Project  Numbers  V32S-R1A,  V32K-13A,  and 
P32K-13A.  Dr.  H.  E.  Scott,  Directorate  of  Technology,  was  the  Air  Force  project  manager, 
The  manuscript  was  submitted  for  publication  on  September  25,  1979. 
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1.0  INTRODUCTION 

Contamination  of  spacecraft  and  satellite  surfaces  and  components  by  exhaust  products 
from  spacecraft  thrusters,  outgassing  products,  and  other  contamination  sources  are  of 
considerable  concern  to  the  spacecraft  designers.  This  concern  is  attributed  to  the 
development  of  more  sophisticated  spacecraft  systems  designed  to  perform  multiple, 
noninterfering  missions  for  longer  periods  of  time  in  space.  Contamination  and  subsequent 
degradation  of  sensitive  optical  components  can  compromise  mission  objectives.  It  is  for 
these  reasons  that  effects  of  condensed  gases  on  optical  surfaces  has  received  widespread 
interest. 

Much  of  the  effort  in  the  past  studies  of  the  radiative  properties  of  cryodeposit-coated 
surfaces  and  materials  has  been  experimental  in  nature.  The  quantity  of  experimental  data 
recently  generated  for  condensed  gases  has  necessarily  emphasized  the  experimental  portion 
of  this  subject.  The  theoretical  and  analytical  side  of  this  study,  although  initiated  earlier  for 
determining  refractive  and  absorption  indices  (n’s  and  k’s)  from  transmission 
measurements,  has  been  lagging.  Further  extension,  therefore,  of  various  mathematical 
models  was  deemed  plausible  for  determining  n’s  and  k’s  of  deposited  films  and  of  optical 
substrate  materials  using  both  transmission  and  reflectance  measurements.  In  addition, 
from  the  transmission  measurements  of  the  pure,  condensed  gases  and  of  mixtures  of  gases, 
efforts  have  been  undertaken  to  try  and  determine  the  mole  fractions  of  various  species  of 
gases  condensed  in  a  matrix.  These  fractions,  determined  by  relative  absorption  band 
strengths,  were  compared  to  mole  fraction  values  determined  experimentally. 

.  The  overall  objectives  of  this  work  are  to  provide  improved  analytical  methods  for 
determining  mole  fractions  of  condensed  rocket  exhaust  gases  and  their  optical  properties 
(n,k)  from  spectroscopic  transmittance  or  reflectance  data.  To  accomplish  these  overall 
objectives,  the  following  specific  areas  were  investigated: 

1 .  Determination  of  the  mole  fractions  of  cryodeposit  constituents. 

2.  Determination  of  the  optical  indices  n  and  k  of  a  medium  from  non-normal 
reflectance  measurements  of  a  single  interface  of  the  medium  with  either  air  or  a 
vacuum. 

3.  Calculation  of  the  absorption  index,  k,  of  a  material  from  known  values  of  the 
refractive  index,  n. 

4.  Obtaining  the  optical  indices  of  a  thin  film  from  transmittance  measurements  of 
the  film  on  a  thick  substrate. 
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5.  Finding  the  optical  indices  of  a  thin  film  from  reflectance  measurements  of  the 
film  on  a  thick  substrate. 

6.  Measurement  of  the  optical  indices  of  substrates  at  cryogenic  temperatures. 

7.  The  use  of  the  molar  refractivity  of  materials  in  mole  fraction  determinations. 

Progress  in  accomplishing  these  items  is  described  in  the  following  sections  of  this  report. 

2.0  MOLE  FRACTIONS  OF  MOLECULAR  SPECIES  IN 
HOMOGENEOUS  CRYODEPOSITS 


2.1  BASIC  THEORY 

Several  spectroscopic  methods  for  finding  concentrations  of  molecular  species  in 
cryodeposits  were  examined  by  one  of  the  authors,  K.  F.  Palmer,  during  the  1976  USAF- 
ASEE  Summer  Faculty  Research  Program  at  AEDC.  The  most  promising  method  is  based 
on  the  assumption  that  the  number  density  and  the  molar  concentration  of  the  i'h  molecular 
species,  c,,  is  directly  proportional  to  the  band  integral  if  species  i  is  infrared  (IR)-active 
(Ref.  I): 


/ kjMdK,  i.e..  c.  cc  fki  (r^)  d V/ 
Band  j  Band  j 


(1) 


Here,  c,  has  the  dimensions  of  moles  per  unit  volume  and  k,(r)  is  the  contribution  of  species  i 
to  the  absorption  index  of  the  homogeneous  medium  at  the  wavenumbers,  v ,  surrounding 
the  j,h  absorption  band  characteristic  of  molecular  species  i.  This  assumption  implies  that 
changes  in  the  intermolecular  forces  should  not  affect  the  integral  of  k,  over  the  absorption 
region  even  though  the  shape  of  the  absorption  index  curve  may  be  altered.  It  is  known  that 
the  formation  and  breaking  of  hydrogen  bonds  can  cause  anomalous  intensity  variations  in 
intramolecular  absorption  bands  (Ref.  2).  Thus,  especially  for  H20  molecules,  Eq.  (1)  may 
need  modification  in  some  instances  (see  Section  2.2). 

A  second  assumption  is  that  tor  k(ri)  of  a  given  homogeneous  material, 

N 

fcM  =  Z 

£=i 
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when  there  are  N  molecular  species  in  the  material  (Ref.  2).  Therefore, 

ni 

=  k(v)  -  £  kj(v)  ^ 

(V. 

where  the  summing  term  does  not  include  the  f=i  term  and  represents  the  effects  of 
"background”  absorption  on  the  absorption  attributed  to  molecular  species  i.  In  practice, 
the  determination  of  the  background  is  somewhat  ambiguous,  and  variations  in  background 
values  can  have  noticeable  effects  upon  the  results  (see  Section  2.2). 


The  molar  concentration  of  any  molecular  species  i,  c,,  in  a  sample  can  be  calculated  if 
the  ratio  R,  =  c/csi  is  known.  The  "standard”  molar  concentration,  cs„  is  the  previously 
measured  concentration  of  species  i  in  a  “standard”  material  that  usually  contains  only  the 
i'h  molecular  species  and  is  at  the  same  temperature  as  the  sample.  For  IR-active  molecular 
species  Eq.  (1)  gives 


R, 


(4) 


where  k,(v)  is  given  by  Eq.  (3)  and  kSI(j>)  is  the  analogous  quantity  for  the  standard  material. 
Note  that  the  band  shapes  in  the  k,(i<)  and  kM(j>)  curves  may  not  be  the  same  and  the  domains 
of  integration  may  not  be  identical. 

Once  the  molar  concentrations  of  the  IR-active  and  inactive  molecular  species  in  a 
material  are  known,  the  mole  fraction  of  species  i,  f„  can  be  calculated  from 

fi  -  Jl  c£  <5> 

It- 1 

where  N  is  the  total  number  of  IR-active  and  inactive  species.  Of  course. 


f.  =  1 


=  i 


(6) 
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and  if  all  but  one  of  the  molar  concentrations  are  known,  the  unknown  molar 
concentration,  say,  fN  where  species  N  is  IR-inactive,  could  be  found  from 

r«  -  '  -  L  f,  <7> 


The  molar  concentrations  of  species  in  materials  are  related  to  the  densities  of  the 
materials.  If  the  standard  material  for  species  i  contains  only  species  i  molecules,  the  density 
of  the  standard  material,  psl,  is 


Ps  i  = 


M. 


(8) 


where  M,  is  the  molecular  mass  of  species  i  with  the  dimensions  of  mass/mole.  If  Eq.  (8)  is 
used  in  the  first  equation  of  Eq.  (4),  the  molar  concentration  of  species  i  in  another  material 
can  be  found  from 


i 


ci  =  P.s.W 


(9) 


However,  R,  can  be  found  from  the  band  integrals  [second  equation  of  Eq,  (4)]  only  for  IR- 
active  species;  R,  for  IR-inactive  species  must  be  obtained  by  other  means  such  as  Raman 
spectroscopy,  or  the  molar  refract ivities  (see  Section  8.0). 

If  only  IR  spectra  are  available,  the  mole  fractions  of  the  molecular  constituents  in  a 
given  material  can  be  found  if  at  most  one  IR-inactive  species  is  present,  whose  identity  is 
known  by  other  methods,  such  as  Raman  spectroscopy.  In  the  case  that  no  IR-inactive 
species  is  in  the  material,  substitution  of  Eq.  (9)  into  Eq.  (5)  gives  the  mole  fraction  of 
species  i  as 


F. 

1 


pst  Vme 


(10) 


The  Appendix  shows  that  for  small  relative  errors  AR|/R,  and  Aes(/gsf  in  Rf  and 
respectively,  the  relative  error  in  f]  is  approximately 
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Af./f.  =  .\R./R.  -  Apv./psl  -  E  h  Aps£/psg) 

E=i 


(id 


Note  in  Eq.  (1 1)  that  if  the  relative  errors  in  the  densities,  qs p,  are  comparable  to  those  of  the 
concentration  ratios,  Rf,  the  influence  of  the  density  errors  and  the  influence  of  the  ratio 
errors  upon  the  calculated  mole  fraction  of  a  species  are  also  comparable. 

If  molecular  species  N  is  the  only  lR-inactive  species  present  in  some  material,  Eq.  (5) 
can  be  written 


N-l 


■  V'(E  H  *  «n) 

e-i 


(12) 


or 


N-l 


fi  =  (ps,R.'M,)/  y  Lp* iVMrcN 


when  Eq.  (9)  is  used  for  the  molar  concentration  of  the  IR-active  species.  The  molar 
concentration  of  the  I R-in active  species  N,  Cn»  cannot  be  found  from  Eq.  (9)  because  Rn  is 
unknown.  However,  RN  may  be  calculated  from  the  average  density,  q,  of  the  sample 
material: 

N  N-l 

p  =  £  cgMg  =  E  eg  Mg  *  ckM>j 

l’=  1  &=l 

N-  l 

Here,  cf  Mf  =  CsfR,'Mf  =  es(Rf,  from  Eq.  (8).  Thus,  cN  =  (e  -  E  QsfRf)/MN  and 
Eq.  (12)  is 


f=  i 


The  relative  error  of  f,  is  approximately 


N-l 

Af/f;  '  ARj-'R;  +  Apxj/psi  *-  -Sp/p  -  £  fg  (Mg'Hlj,  -  lXARg'Rg  +  Ap^/p^g  -  Ap/p>  <]4) 

l=\ 


when  the  errors  in  Rf,  and  q  are  small  (see  the  Appendix).  Again,  the  relative  errors  of 
the  densities  can  be  just  as  important  as  the  relative  errors  of  the  concentration  ratios  in 
causing  errors  in  the  calculated  mole  fractions. 
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In  summary,  the  mole  fractions  of  molecular  constituents  in  a  sample  material  can  be 
found  from  infrared  spectra  using  either  Eq.  (10)  or  (13)  if  the  densities  of  the  sample  and 
the  standard  materials  are  known  and  if  the  ratio  R  of  the  sample  molar  concentration  to  the 
molar  concentration  of  the  standard  material  are  measured  for  all,  or  all  but  one, 
constituent  molecular  species.  For  lR-active  species  the  concentration  ratios  R  are  assumed 
identical  to  ratios  of  the  appropriate  band  integrals  (Ref.  1).  Also,  accurate  determination  of 
the  mole  fractions  of  species  depends  as  much  on  the  accuracy  of  the  density  values  as  on  the 
accuracy  of  the  band  integral  ratio  values. 

2.2  EXPERIMENTAL  RESULTS 

Table  1  contains  the  results  of  the  application  of  Eqs.  (10)  and  (13)  to  AEDC  cryodeposit 
data  previously  analyzed  by  W.  B,  Newbolt  during  the  1977  USAF-ASEE  Summer  Faculty 
Research  Program  at  AEDC.  In  the  two  columns  next  to  the  last  one  the  mole  fractions 
calculated  in  the  present  analysis  are  compared  to  the  mole  fractions  found  from  chemical 
analysis  of  the  residual  gas  mixture  used  to  form  the  cryodeposits.  Note  that  the  calculated 
mole  fractions  are  not  identical  to  the  ratios  of  the  band  integrals.  In  the  last  column  the 
mole  fractions  calculated  by  Eq.  (10)  or  (13)  from  the  band  integral  ratios  of  Newbolt  are 
given.  In  most  cases  Newbolt’s  choice  of  band  integrals  differed  from  those  in  the  table,  and 
there  were  no  allowances  made  for  any  background  absorption. 

The  agreement  of  the  present  mole  fraction  values  with  those  from  chemical  analysis  is 
within  15  percent  for  the  two-component  cryodeposits,  but  is  much  worse  for  the  Plume  I 
film  composed  of  a  simulated  plume  mixture  of  gases.  In  particular,  the  concentration  of 
water  (H20)  molecules  was  calculated  to  be  much  less  than  was  found  chemically.  This 
discrepancy  is  most  likely  due  to  the  anomalous  behavior  of  H20  absorptions  caused  by 
hydrogen  bonding,  an  effect  that  depends,  in  part,  on  the  concentration  of  H20  molecules. 
For  pure  H20,  the  absorption  band  intensities  of  intramolecular  vibrations  of  condensed 
phases  are  approximately  an  order  of  magnitude  greater  than  the  corresponding  vapor  phase 
intensities  attributable  mainly  to  the  increased  hydrogen  bonding  in  the  condensed  phases 
(Ref.  2).  Although  more  study  on  the  effects  of  hydrogen  bonding  in  cryodeposits  is 
necessary,  it  should  be  noted  that  the  band  integral  of  the  librational  band  of  liquid  water, 
Pi,t>,  is  very  nearly  temperature  independent  and  so,  presumably,  is  nearly  independent  of  the 
local  molecular  environment  (Ref.  3).  It  might  be  possible  to  use  the  librational  band  of  H20 
in  cryodeposits  as  a  more  accurate  index  of  H20  concentration.  To  do  this  the  domain  of  the 
IR  data  must  be  extended  down  to  about  500  cm-1;  this  can  only  be  done  by  obtaining 
accurate  measurements  of  the  optical  properties  of  the  substrate  at  cryogenic  temperatures 
(see  Section  7.0). 
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In  the  calculations,  a  consistent  set  of  cryodeposit  densities  is  used  throughout.  In  all 
cases,  however,  better  agreement  with  the  mole  fraction  determined  from  chemical  analysis 
is  possible  if  slightly  different  values  of  the  densities  are  used  in  each  separate  calculation. 
The  densities  given  in  Table  1  for  the  dinitrogen/water  mixture  (N2/H20),  carbon  monoxide 
(CO),  and  Plume  I  cryodeposits  as  well  as  that  of  pure  N2  (which  was  not  used  in  the 
calculations)  are  those  found  during  the  original  transmission  measurements.  The  density  of 
pure  carbon  dioxide  (CO2)  is  the  value  determined  by  Arnold  et  al.  (Ref.  4),  that  of  H20  is 
the  mean  value  of  several  H20  cryodeposits  measured  at  AEDC,  and  the  densities  of  the 
remaining  mixtures  represent  the  mean  of  the  tabulated  densities  of  the  standard  (pure) 
cryodeposit  of  the  constituent  species,  weighted  according  to  their  approximate  molar 
concentrations. 

Except  for  those  marked  by  asterisks,  the  band  strengths  (integrals)  given  in  Table  1  are 
found  by  numerical  integration  of  the  data  over  the  indicated  wavenumber  domains.  If  the  k 
value  at  either  endpoint  of  the  domain  was  not  zero,  a  background  absorption  was 
computed  and  subtracted  from  the  integrals  to  give  the  tabulated  band  strength,  S. 

There  is  a  degree  of  arbitrariness  in  choosing  the  limits  of  integration  and  accounting  for 
extraneous  absorptions.  A  “moment”  method  for  finding  the  baseline  and  the  limits  of 
integration  of  unblended  bands  has  been  shown  by  Jiang  et  al.  (Ref.  5)  to  reduce  the 
experimental  uncertainties  in  the  computation  of  band  strengths,  and  it  is  recommended  that 
the  “moment”  method  be  tried  and  evaluated  in  future  band  strength  calculations.  Blended 
absorptions  were  decomposed  by  a  nonlinear  least-squares  fit  of  the  data  to  two  or  more 
Lorentzian  line  shapes.  For  instance,  the  v2  band  of  H20  is  badly  blended  by  lattice 
vibration  absorptions  in  the  pure  H;0  cryodeposit  and  is  blended  by  other  absorptions  in  the 
Plume  I  cryodeposit.  Lattice  bands  also  contribute  to  the  absorption  of  the  1—0  band  of 
12CO  and,  to  a  lesser  extent,  the  vj  band  of  l2CC>2  (Ref.  5).  For  these  reasons,  absorptions 
attributable  to  lJCO  and  l?C02  were  used  whenever  possible,  although  the  much  lower  band 
intensities  of  the  nC  isotopic  species  are  subject  to  larger  experimental  errors.  For  the  Plume 
1  cryodeposit,  only  the  ^3  band  of  l2CO  was  sufficiently  “clean”  of  extraneous  absorptions 
so  that  its  band  strength  could  be  found  by  numerical  integration  of  the  original  data.  The 
other  band  strengths  in  the  Plume  I  cryodeposit,  and  those  of  the  standard  cryodeposits,  are 
computed  from  Lorentzian  fits  because  of  the  blends  cited  above.  It  should  be  noted  that  the 
Lorentzian  band  strength  for  the  1  —0  band  of  l2CO  is  greater  than  the  band  strength  found 
by  direct  integration,  because  of  the  large  “wings”  of  the  Lorentzian  line  shape.  This 
suggests  that  a  Gaussian  or  Voight  profile  might  better  approximate  the  line  shapes  of  the 
data. 
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3.0  OPTICAL  INDICES  OF  A  MEDIUM  FROM  NON-NORMAL 
REFLECTANCE  MEASUREMENTS  OF  A  SINGLE  INTERFACE 

The  optical  constants  n(v)  and  k(v)  al  wavenumber  i>  of  a  material  can  be  obtained  from 
reflection  measurements  of  a  light  beam  at  the  interface  of  the  material  with  a  vacuum  or 
air.  Kramers-Kronig  (KK)  techniques  allow  both  n  and  k  to  be  found  from  the  reflectance  of 
a  beam  of  electromagnetic  radiation  at  the  interface  using  measurements  at  just  one  angle  of 
incidence.  Of  course  the  material  must  have  dimensions  and  a  geometry  that  will  not  allow 
extraneous  radiation  to  enter  the  spectrometer;  precautions  are  especially  necessary  at 
wavenumbers  where  the  material  is  nearly  transparent.  For  non-normal  incidence  the 
components  of  the  electric  and  magnetic  fields  parallel  and  perpendicular  to  the  plane  of 
incidence  have  reflectances  that  become  measurably  different  at  beam  incident  angles 
greater  than  about  10  deg  (Ref.  6). 

The  complex  reflectivity  if  r(v)  =  [R(v)],/;  exp  [i0(jO],  where  the  circumflex  denotes 
complex  quantities,  i  is  (-1)1/2,  R  is  the  measured  reflectance,  and  $  is  the  phase  shift  of  the 
beam  upon  reflection.  Two  KK  dispersion  relations  between  R  and  <j>  can  be  written  in  which 
one  of  the  quantities  is  evaluated  at  a  given  wavenumber  in  terms  of  an  integral  of  the  other 
quantity  over  all  wavenumber  space.  In  our  case  the  reflectance,  R,  is  measured 
experimentally  and  the  phase  shift,  </>,  can  be  found  from 


fn  [R(„  h]V2  fh,- 


2  *  2 
v  —  v 


(15) 


where  P  denotes  the  Cauchy  principle  value.  The  optical  constants  are  obtained  from 

n  = 

and 

k  -  SQ/n 

where 

Q  =  (1  -  R)  cos  6  /  (1  +  H  -  2R*  cos  <f>) 

and 

S  =  2  R  sin  cos  8  }  (l  +  R  —  2  II ^  cos  cO 

for  a  beam  angle  of  incidence  &.  Note,  however,  that  for  non-normal  incident  angles,  the 
reflectances  of  different  polarization  components  of  the  beam  will  vary.  It  is  important  that 
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the  reflectance  of  only  one  polarization  component  be  used  in  these  equations  when  6  >  > 
10  deg. 

The  finite  range  of  measured  reflectance  data  means  that  extrapolations  of  R-values  to 
v  =  0  and  to  v  =  co  are  necessary.  In  the  present  case  the  R-values  are  assumed  to  be 
constant  outside  the  measured  data  domain  with  R(v>  =  R(t<L)  for  0  <  v  <  yL,  and  R(j>)  = 
R(i^h)  for  v  >  vy\. 

The  errors  introduced  by  various  extrapolations  of  reflectance  data  and  other 
experimental  uncertainties  have  been  discussed  by  Bowlden  and  Wilmshurst  (Ref.  7)  and 
Hale  et  al.  (Ref.  8).  The  former  authors  have  studied  the  effects  of  neglecting  all 
contributions  to  the  integral  of  Eq.  (5)  outside  the  domain  of  the  data.  This  is  equivalent  to 
assuming  that  R  is  1  when  0  <  v  <  pL  or  v  >  Using  comparisons  to  theoretical 
dispersion  curves,  they  show  that  the  k-values  from  the  KK  analysis  of  an  abosrption  band 
are  within  a  few  percent  of  the  theoretical  values  as  long  as  the  data  limits  are  sufficiently  far 
from  the  band  center,  say,  greater  than  10  band  half-widths.  When  they  varied  the  limits  of 
integration,  which  is  analogous  to  changing  the  extrapolation  conditions  of  actual 
reflectance  data,  the  most  noticeable  effect  they  found  was  some  (asymmetric)  distortion  of 
the  band  shape,  a  result  in  accord  with  the  authors’  experience  with  liquids.  This  effect  was 
minimized  for  a  given  band  centered  at  p0  when  the  limits  of  integration  were  chosen  so  that 
p0/pl  *»  ph/vo-  Under  these  conditions  the  errors  caused  by  neglecting  the  high  and  low 
wavenumber  regions  cancel  almost  completely,  except  for  wavenumbers  within  a  small 
region  about  the  band  center.  There  the  errors  are  additive  but  very  small.  Thus,  it  appears 
that  the  k  values  determined  by  the  KK  method  are  relatively  insensitive  to  the  details  of 
extrapolations  used  outside  the  data  domain,  and  it  is  the  authors’  experience  that  the  n 
values  determined  by  the  KK  method  are  even  less  sensitive. 

Hale  et  al.  (Ref.  8)  have  shown  that  when  the  reflectance  is  analyzed  by  KK  techniques, 
the  effect  of  random  relative  errors  in  the  IR  reflectance  values  of  liquid  water  cause 
approximately  the  same  values  of  relative  error  in  all  n  values  and  in  the  k  values  within  an 
absorption.  The  random  relative  errors  for  k  values  become  large  at  wavenumbers  far  from 
absorptions.  These  results  should  apply  to  work  on  other  dielectrics  because  of  similar 
reflectance  and  error  values. 

The  influence  of  systematic  errors  in  the  reflectance  upon  the  optical  constants  from  KK 
analysis  is  felt  to  be  minimal.  Consider  the  special  case  in  which  all  “true”  reflectance 
values,  R,  are  changed  by  some  systematic  error  to  R'  =  /3R,  where  &  is  a  constant.  Now, 
Eq,  (15)  when  integrated  by  parts  is 


fn^fl(  o  ) 
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[f  R'  were  used  in  the  integral  instead  of  R,  the  differential  would  become  d(lnvTT)  = 
d(ln\/£fl£)  =  d(lnV^)  +  d(lnVR)  =  d(lnVK)  are  the  phase  shift,  <t> ,  would  remain  unchanged 
from  the  “true”  value.  In  other  instances  of  systematic  error  where  0  is  a  slowly  varying 
function  of  wavenumber  compared  to  In^R  the  changes  in  4>  would  be  very  small. 

Ahrenkiel  (Ref.  9)  has  shown  that  subtractive  Kramers-Kronig  (SKK)  dispersion 
relations  reduce  the  errors  in  the  optical  constants  associated  with  extrapolations.  Essentially 
SKK  values  of  the  optical  constants  are  made  relative  to  a  known  value  of  that  constant  at  a 
wavenumber  v0)  which  is  within  the  wavenumber  domain  of  the  data.  For  reflectance 
measurements  from  a  single  air-material  or  vacuum-material  interface,  the  appropriate  SKK 
relation  for  the  phase  shift  is 


©(v) 


V 


l’n  It  (v  0  civ  ' 


(16) 


Here,  R  is  the  reflectance  of  a  single  component  of  polarization  that  is  usually  the 
component  perpendicular  to  the  plane  of  incidence  if  the  angle  of  incidence,  6,  is  not  zero. 
The  optical  constants  can  be  found  from  the  above  relations. 


Computer  programs  “KKANR”  and  “SKKR”  (written  by  K.  F.  Palmer)  incorporate 
Eqs.  (15)  and  (16),  respectively,  for  KK  and  SKK  analysis  of  reflectance  data.  The  input  data 
to  the  programs  can  be  arbitrarily  spaced  according  to  wavenumber.  The  integration  is  done 
numerically  by  Simpson’s  rule  except  near  the  singular  points  and  where  three  consecutive 
data  points  are  not  equally  spaced  in  wavenumber.  In  these  cases  the  numerator  of  the 
integrand  is  approximated  by  a  quadratic  function  Fitted  to  three  consecutive  data  points, 
and  the  integration  is  computed  analytically. 

Both  “KKANR”  and  “SKKR”  have  been  used  to  analyze  the  ice  reflectance  data  of 
Schaaf  (Ref.  10).  The  agreement  between  the  two 'programs  for  the  optical  constants  is 
excellent  if  some  care  is  exercised  is  choosing  for  the  latter  program,  such  as  using  a  v0 
value  that  is  far  from  any  absorptions.  Also,  the  values  determined  by  these  programs  for  all 
n  values  and  those  k  values  surrounding  molecular  absorptions  are  within  one  percent  of  the 
values  given  by  the  KK  analysis  of  Schaaf.  However,  the  k  values  far  from  absorptions  are 
sometimes  substantially  different  from  those  of  Schaaf  and  may  even  be  negative.  Part  of 
the  discrepancy  may  be  due  to  the  difference  in  algorithms  for  obtaining  the  contribution  to 
the  integrals  about  a  singular  point.  The  authors  feel  the  present  algorithm  is  more  accurate, 
however,  and  suggest  further  investigation  of  the  matter. 
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Kramers-Kronig  dispersion  relations  can  be  written  for  the  real  and  imaginary 
components,  or  the  modulus  and  phase  angle,  of  many  well-behaved  complex  functions.  In 
particular,  the  KK  equations  for  n  and  k  (the  real  and  imaginary  components  of  the  complex 
refractive  index  n)  can  be  used  to  find  k  if  n  is  known  over  a  wide  wavenumber  region,  and 
vice-versa.  The  relations  can  also  be  used  to  check  the  internal  consistency  of  optical  data. 

The  KK  relations  for  n  and  k  are 


n(i') 


v  1  ft.- ')  dt/ ' 
*> 

v  "  —  v  “ 


07) 


and 


kit') 


n  dv  ' 

0  *1 
v  —  V 


(18) 


where  n(oo)  is  the  contribution  to  the  index  of  refraction  as  v  goes  to  infinity  (Ref.  4).  The 
corresponding  SKK  relations  are  easily  shown  to  be 


and 


n(i/)  =  nh'0) 


_ i* V  I|/  h  <U/ 

(|/^  —  v  '*)  (t/“  -  v'~) 


(19) 


k(i/i  = 


k(v  )  2C'2  -  r/2l 

- 2_  +  ...  .  n.  .  .  — 

V  !7 


ntt'  h  dv' 


‘•v“  —  v  '■")  —  V  ^ 


(20) 


where  n(r0)  and  k(v0)  are  found  at  v  =  v0  by  a  separate  measurement.  A  computer  program 
utilizing  Eq.  (19)  is  presently  in  use  at  AEDC  (Ref.  1 1).  Computer  programs  “KKN2K”  and 
“SKKN2K”  have  been  written  (by  K.  F.  Palmer)  to  calculate  k  from  known  n  values  using 
Eqs.  (18)  and  (20),  respectively.  Outside  the  data  domain  i>|_  <  v  <  t'n  the  values  of  k  are 
assumed  to  be  k(vH)  for  the  v  >  and  k(iq )  for  0  <  v  <  vL. 
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Program  “SKKK2N”  has  also  been  written  (by  K.  F.  Palmer)  to  obtain  n  values  from  k 
values.  As  a  check  on  the  internal  consistency  of  programs  “SKKK2N”  and  “SKKN2K”, 
they  were  run  in  sequence.  The  AEDC-observed  values  of  k(*)  in  the  range  from  700  cm  - 1 
to  3,700  cnr1  (including  the  finely  spaced  points  in  the  absorption  bands)  for  CO2  were 
used  in  SKKK2N  to  obtain  the  n(»<)  spectrum.  These  n(v)  values  were  then  used  as  input  for 
the  SKKN2K  program,  yielding  a  k(»»)  spectrum  (k-calculated)  that  could  then  be  compared 
with  the  original  AEDC  values.  The  differences  between  the  original  k-observed  and 
k-calculated  values  were  substantial— consistently  less  than  0.2  except  at  the  peak  of  the 
band  of  C02,  where  the  difference  was  approximately  0.3,  and  near  the  lower  end  of  the 
observed  spectrum  (710  to  1030  cnr')  where  the  differences  were  as  much  as  0.565. 

There  are  several  plausible  explanations  for  these  discrepancies,  including  coding  errors 
in  the  computer  programs,  numerical/analytical  errors  (such  as  formula  derivations),  and 
inaccurate  data.  It  is  encouraging  that  the  error,  (k-observed)  -  (k-calculated),  appears  to  be 
quite  regular.  Indeed,  the  values  of  k-observed  and  k-calculated,  as  well  as  their  difference, 
were  plotted  and  compared.  The  V3  band  of  C02  in  particular  had  the  same  shape  and 
location,  but  appeared  to  be  shifted  downward  by  approximately  0.2.  This  leads  us  to 
believe  that  there  may  be  a  programming  or  numerical/analytical  error,  and  thus  further 
investigation  of  the  consistency  of  the  SKK  relations,  and  the  corresponding  programs,  will 
begin  here. 


5.0  OPTICAL  CONSTANTS  OF  A  THIN  FILM  ON  A  THICK 
SUBSTRATE  FROM  TRANSMITTANCE  MEASUREMENTS 

At  present  at  AEDC,  the  optical  constants  of  thin  films  deposited  on  thick  substrates  are 
determined  by  a  separate  calculation  at  each  wavenumber  using  the  dependence  of  the 
transmittance  of  a  beam  through  the  film  and  substrate  upon  the  thickness  of  the  film.  A 
lamellar  model  of  the  film  and  substrate  (Ref.  12),  which  includes  interference  effects,  is 
treated  by  a  nonlinear  least-squares  procedure  to  extract  the  optical  constants  at  a  given 
wavenumber.  One  of  the  underlying  assumptions  of  this  method  is  that  the  optical 
properties  of  the  film  do  not  change  with  thickness. 

Meada  et  al.  (Ref.  13)  have  shown  that  KK  relationships  can  be  written  for  the  modulus 
and  phase  angle  of  the  complex  transmission  coefficient,  t,  of  a  single  film  thickness  on  a 
substrate,  where  t(?)  =  T(y)l/2  e1^  and  T(r)  is  the  measured  transmittance  at  wavenumber 
v.  The  KK  relation  of  interest  is 

00 

*(!/)  =  ■£  P  f  -fri-W  V  +  2r,ud  (21) 
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where  d  is  the  film  thickness.  The  corresponding  SKK  equation  is 


4>  (v^ 


(!n  T  (i/J)  dv' 


lu2-w'2)  Wl~v'2) 


(22) 


<£(r0)  being  the  phase  angle  at  v0  determined  from  a  separate  measurement  at  that  film 
thickness.  Either  equation  can  be  used  to  evaluate  o  of  a  single  thickness  of  film,  and  both 
optical  constants  n  and  k  can  be  found  for  that  thickness  by  a  nonlinear  interpolation  of  the 
complicated  functions  relating  n,  k,  T,  and  4>  (Ref.  13). 


A  computer  program,  which  analyzes  transmittance  data  at  a  single  film  thickness 
according  to  Eq.  (22),  was  applied  to  AEDC  data  of  a  CO  cryodeposit  with  a  film  thickness 
near  the  midrange  of  the  observed  thicknesses.  In  most  instances  this  SKK  program  was 
successful  in  converging  to  the  same  n  and  k  values  found  by  using  the  dependence  of  the 
transmittance  upon  film  thickness.  Sometimes,  however,  especially  near  molecular 
absorptions,  more  than  one  physically  reasonable  pair  of  n  and  k  values  will  give  identical 
values  of  T  and  <f>,  and  the  two  methods  will  not  agree.  At  present  potential  algorithms  for 
finding  the  “best”  values  of  n  and  k  are  being  tested.  It  may,  for  example,  be  necessary  to 
analyze  data  from  two  (or  more)  film  thicknesses  simultaneously,  so  that  by  comparison  of 
the  possible  solutions  for  both  thicknesses,  the  unwanted  solutions  might  be  eliminated. 

6.0  OPTICAL  CONSTANTS  FROM  REFLECTANCE  MEASUREMENTS 
OF  A  THIN  FILM  ON  A  SUBSTRATE 


Radiation  reflected  from  a  thin-film  cryodeposit  on  a  metal  substrate  has  been  analyzed 
by  using  the  dependence  of  the  reflectance  upon  film  thickness  at  a  given  wavenumber  (Ref. 
4).  As  in  the  analysis  of  transmittance  data,  a  lamellar  model  of  the  film  and  the  substrate 
was  used,  and  it  was  assumed  that  the  optical  constants  of  the  film  are  independent  of  film 
thickness. 

From  the  results  given  in  Section  5.0,  it  seems  feasible  that  usable  KK  relations  exist  for 
the  modulus  and  phase  angle  of  the  complex  reflection  coefficient  r(i>)  =  R(v)l/2  e"^*, 
where  R(p)  is  the  measured  absolute  reflectance  (or  the  reflectance  relative  to  a  bare 
substrate)  of  a  single  film  thickness.  The  main  criterion  for  determining  the  usefulness  of  a 
KK  relation  is  that  the  value  of  the  desired  optical  property  be  finite  and  determinate  as  the 
wavenumber  approaches  zero  or  infinity.  At  present  the  value  of  0(0°)  of'r'is  not  known, 
although  a  preliminary  investigation  indicates  the  value  is  very  likely  zero.  If  a  suitable  0{<») 
value  is  found  in  future  work,  computer  programs  similar  to  the  one  being  tested  (Section 
5.0)  can  result,  giving  a  KK  method  of  obtaining  the  optical  constants  n  and  k. 
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7.0  OPTICAL  CONSTANTS  OF  SUBSTRATE  MATERIALS 
AT  CRYOGENIC  TEMPERATURES 

In  the  lamellar  model  of  a  thin  film  on  a  substrate,  the  optical  constants  of  the  film  can 
be  found  if  those  of  the  substrate  are  known  (see  Sections  5.0  and  6.0).  There  are,  however, 
few  reliable  data  on  the  cryogenic  optical  properties  of  substrate  materials  in  the  literature, 
and  “in-house”  measurements  must  be  made.  At  present,  measurements  of  the  optical 
constants  of  germanium  at  20  and  80  K  have  been  made  in  the  range  from  700  to  3,700  cm  - 1 
(Ref.  12).  Section  2.0  describes  the  need  to  observe  molecular  absorptions  of  cryodeposits 
beyond  the  high  and  low  limits  of  this  wavenumber  region.  This  can  be  accomplished  if  the 
range  of  germanium  measurements  is  expanded  and/or  data  are  obtained  for  other  substrate 
materials. 


The  analysis  of  the  transmittance  or  reflectance  of  two  substrate  thicknesses  is  the  most 
direct  w'ay  of  determining  the  optical  properties  of  the  substrate,  but  this  analysis  may  be 
expensive  or  otherwise  inconvenient.  A  KK  analysis  of  a  single  thick  film  of  substrate 
material  is  not  possible  because  of  the  loss  of  phase  information.  Another  alternative  for 
finding  substrate  optical  properties  is  the  analysis  of  reflectance  measurements  from  a  single 
interface  of  substrate  material  through  an  appropriate  KK,  or  SKK,  relation  (see  Section 
3.0).  However,  shaping  the  sample  into  a  geometry  such  as  a  light  horn  (Ref.  14),  that  is, 
suitable  for  the  elimination  of  stray  reflections,  may  not  be  practical  for  some  materials, 
especially  for  transparent  wavenumber  regions. 

8.0  MOLAR  REFRACTIVITY  OF  MATERIALS  -  DETERMINATION 
OF  MOLECULAR  SPECIES  CONCENTRATIONS  IN  MATERIALS 
CONTAINING  MORE  THAN  ONE  INFRARED-INACTIVE  SPECIES 


The  molar  refractivity,  A,  of  a  material  containing  a  single  nonpolar  molecular  species 
is,  at  wavenumbers  far  from  absorptions. 


n  +  2 


M_  n2  -  [ 
p  n2  +  2 


(23) 


where  NA  is  Avogadro’s  number,  «  is  the  molecular  polarizability,  N  is  the  number  density, 
M  is  the  molecular  mass  per  mole,  and  e  is  the  mass  density  of  the  material  (Ref.  6).  If  there 
are  N  chemically  inactive  molecular  species  in  a  homogeneous  mixture,  with  density  g  and 
refractive  index  n,  the  molar  refractivity  of  the  mixture  is  approximately 

N  N 

A  =  I  c.At/  l  c. 

i=l  i=l 
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where  c,  is  the  molar  concentration  (moles  per  unit  volume)  of  species  i  in  the  mixture  and  A, 
is  its  molar  refractivity.  Equation  (23)  and  the  above  equation  can  be  combined  to  show 


N 


z 


c  M 
9 


n“  ~  1 
n2  +  2 


N 


Z 

i=  l 


r.  M 

i  i 

P  . 


s  J 


(24) 


where  M,  is  the  molecular  mass  per  mole  of  species  i,  gsl  is  the  known  density,  and  nsl  is  the 
known  refractive  index  of  the  “standard”  material  for  species  i. 

In  principle,  the  concentrations  c,  could  be  determined  from  the  values  of  the  densities 
and  indices  of  refraction  measured  at  least  N  different  wavenumbers,  allowing  the 
calculation  of  the  mole  fractions  by  Eq.  (5)  of  Section  2.0.  However,  for  wavenumbers  far 
from  absorptions  the  indices  ordinarily  do  not  change  noticeably  with  wavenumber,  and  a 
numerical  solution  of  Eq.  (24)  may  be  difficult,  or  impossible,  to  achieve  for  a  material  with 
many  molecular  species.  It  appears  feasible,  though,  that  an  interalive  method  incorporating 
Eqs.  (4)  and  (24)  for  the  IR-active  and  inactive  species,  respectively,  could  be  devised  so  that 
the  molecular  concentrations  can  be  obtained  for  materials  containing  more  than  one  IR- 
inactive  molecular  species.  Of  course,  as  for  the  materials  containing  just  one  IR-inactive 
species  (Section  2.0),  the  identities  of  the  inactive  species  must  be  found  from  techniques 
other  than  IR  spectroscopy,  e.g.,  Raman  spectroscopy.  A  matter  of  further  study  will  be  to 
determine  the  effect  of  applying  Eq.  (24)  to  materials  containing  strongly  polar  molecules 
such  as  H2O. 


9.0  SUMMARY 

The  overall  objective  of  this  work  was  to  provide  improved  analytical  methods  for 
determining  mole  fractions  of  condensed  rocket  exhaust  gases  and  their  optical  properties 
from  spectroscopic  transmittance  and  reflectance  data.  Seven  specific  areas  (see  Section  1.0) 
were  investigated  to  further  progress  toward  the  overall  goal,  and  the  results  may  be 
summarized  as  follows. 

1.  The  investigation  made  for  determining  molar  concentrations  indicated  that 
mole  fractions  can  be  found  from  infrared  spectra  if  the  densities  of  the  sample 
and  standard  materials  are  known  and  if  the  ratio  R  of  the  sample  molar 
concentration  to  the  molar  concentration  of  the  standard  material  is  measured 
for  all,  or  all  but  one,  constituent  molecular  species.  It  was  also  shown  that 
accurate  determination  of  the  mole  fractions  of  species  depends  as  much  on  the 
accuracy  of  the  density  values  as  on  the  accuracy  of  the  absorption  band  integral 


\ 
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ratio  values.  Good  agreement  (within  1 5  percent  or  less)  was  obtained  between 
experimentally  determined  mole  fractions  and  the  integrated  band  model  for  a 
mixture  of  two  components.  For  a  four-component  mixture  fair  agreement  was 
obtained  for  the  CO  and  C02  components  (see  Table  1),  whereas  the 
concentrations  determined  for  H20  and  N2  were  not  close. 

2.  The  optical  constants  n(*0  and  k(r)  of  a  material  can  be  obtained  from 
reflectance  measurements  at  the  interface  of  a  material  with  vacuum  or  air.  Two 
types  of  Kramers-Kronig  programs  for  these  calculations  were  developed. 

3.  The  Kramers-Kronig  analysis  for  finding  k(i>)  from  n(y)  values  was  attempted. 
Theoretically,  this  should  have  been  straightforward,  but  poor  results  were 
obtained  using  the  experimental  data  for  C02. 

4.  A  model  was  developed  and  a  program  written  for  determining  the  n’s  and  k’s  of 
a  thin  film  on  a  thick  substrate  from  transmittance  measurements.  In  the  past, 
many  film  thicknesses  were  required  to  achieve  these  results.  This  single 
thickness  model  gave  good  results  for  a  CO  condensed  film  except  in  the  region 
of  absorption  bands  where  more  work  needs  to  be  done. 

5.  A  similar  study  to  the  previous  one  was  carried  out  for  determining  n’s  and  k’s 
from  a  single  thin  film  from  reflectance  measurements.  However,  the  phase 
angle  near  infinity  <£(00)  is  not  known  so  this  model  cannot  be  used  until  this 
value  is  determined. 

6.  The  problem  of  finding  the  optical  constants  of  optical  substrate  materials  at 
cryogenic  temperatures  was  also  investigated  briefly.  It  appears  that  analysis  of 
the  transmittance  or  reflectance  of  two  substrate  thicknesses  is  the  most  direct 
way  of  determining  these  properties. 

7.  The  last  area  investigated  was  the  determination  of  molecular  species 
concentrations  in  materials  containing  more  than  one  infrared-active  species 
from  material  molar  refractivities.  The  final  conclusion  of  this  study  was  that  it 
appears  feasible  that  an  iterative  method  could  be  devised  such  that  the 
molecular  concentrations  could  be  obtained  for  materials  containing  more  than 
one  IR-active  specie. 
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Table  1.  Calculated  Mole  Fractions  of  Molecular  Species  in  Selected  20  K  Cryodeposits 
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APPENDIX 

PROPAGATION  OF  ERRORS  IN  MOLE  FRACTION  DETERMINATIONS 


In  the  following  discussion,  the  only  errors  are  assumed  to  be  in  the  ratios  of  band 
strengths,  R,  and  in  the  densities,  g,  of  the  materials. 


Case  I.  If  there  are  no  IR-inactive  molecular  species  in  the  material  and  Eq.  (10), 


is  used,  then 


Now,  if  i  * 


Also, 


because  of  the  symmetry  of  the  roles  of  qs(  and  Rf  in  Eq.  (10). 
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by  symmetry.  Upon  substitution  of  these  partial  derivatives  into  the  expression  for  df,  one 
obtains,  after  some  algebraic  manipulation, 
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For  small  errors  in  Rrand  qs(,  Eq.  (11)  follows. 


Case  2.  If  only  one  IR-inactive  molecular  species  (labeled  N)  is  in  the  material  and  Eq  (13) 
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by  symmetry.  If  t  =  i, 
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by  symmetry.  Also, 


where 
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After  substituting  these  partial  derivatives  into  the  expression  for  df,  and  doing  some 
algebraic  manipulation,  one  finds 
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For  small  errors  in  Rf,  qs,\  and  e,  Eq.  (14)  follows. 
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NOMENCLATURE 

A  Molar  refractivity 

c  Molar  concentration  of  iIh  molecular  specie 

i 

f  Mole  fraction 

k  Absorption  index 

M  Molecular  mass 

N  Number  density 

Na  Avrogadro’s  number 

n  Index  of  refraction 

P  Cauchy  Principle  value 

R  Measured  reflectance 

Complex  reflectivity 
S  Molecular  band  strength 

T  Complex  transmission  coefficient 

^  Measured  transmittance 

a  Molecular  polarizability 

$  Incidence  angle 

q  Density,  gm/m3 

</*  Phase  shift 

SUBSCRIPTS 

H  High 

i  Species 

L  low 

lib  Librational 
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f  Material 

N  Species  N 

s  Standard 
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